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Abstract 

Résumé. Même si elle est non-intégrable, une équation diflFérentielle peut 
néanmoins admettre des solutions particulières globalement analytiques. Sur 
l'exemple du système dynamique de Kuramoto et Sivashinsky, génériquement 
chaotique et d'un grand intérêt physique, nous passons en revue diverses métho- 
des, toutes fondées sur la structure des singularités, permettant de caractériser 
la solution analytique qui dépend du plus grand nombre possible de constantes 
d'intégration. 

English abstract. Even if it is nonintegrable, a differential équation may 
nevertheless admit particular solutions which are globally analytic. On the 
example of the dynamical System of Kuramoto and Sivashinsky, which is gener- 
ically chaotic and présents a high physical interest, we review various methods, 
ail based on the structure of singularities, allowing us to characterize the an- 
alytic solution which dépends on the largest possible number of constants of 
intégration. 

Mots clefs : ondes solitaires, équations de Briot et Bouquet, Nevanlinna, algorithme 
LLL, équation de Kuramoto et Sivashinsky, équation de Ginzburg et Landau complexe 
cubique unidimensionnelle, genre, troncature. 
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|6 Ver s une approche numé rique de la solution elliptiauel 



6.3 Détermination des degrés de la sous-éauatiorj 



I7 Approches analvtiauesl 

7.1 Équation de Monge- Ampère éauivalentel 



7.3 Construction de la sous-éauation par la série de Laurent! 



1 Situation du problème 



D'après un résultat classique |28j . un système différentiel ne peut présenter un com- 
portement chaotique que si son ordre égale au moins trois. D'autre part, malgré 
l'opposition intuitive entre chaos et analyticité, il n'est pas interdit à un système 
génériquement chaotique de posséder des solutions particulières globalement analy- 
tiques, la seule impossibilité étant que le nombre de constantes d'intégration ne peut 
égaler l'ordre. Le but de cet article est de passer en revue les diverses méthodes 
connues permettant de caractériser (idéalement d'exhiber) cette solution analytique, 
forcément particulière, qui dépend du nombre maximal possible de constantes d'inté- 
gration. Convenons de l'appeler solution analytique générale, par opposition à la 
solution générale dont l'expression analytique n'existe pas. 

L'exemple physique optimal pour illustrer ces méthodes est à notre avis l'équation 
de Ginzburg et Landau cubique unidimcnsionnellc (CGL3) 

iipt +pipxx + q\ip\'^i^ - ijA^ 0, pq^ ^ 0, e C, 7 G 7^, (1) 

pour les raisons suivantes, 

1. elle est génériquement chaotique |21| . 

2. elle admet une limite intégrable — ^{q) — j — (équation de Schrôdinger 
nonlinéaire (NLS)) dont tous les éléments d'intégrabilité sont connus analytique- 
ment, 

3. les phénomènes physiques qu'elle modélise sont très importants : propagation du 
signal dans une fibre optique ^ , intermittence spatio-temporelle en turbulence, 
supraconductivité, etc. 

Nous nous contenterons d'étudier ici une équation plus simple, celle de Kuramoto et 
Sivashinsky (KS), déduite de celle pour la phase ip ~ arg?/' du champ ip de CGL3 en 
prenant une certaine limite (25ii if 9j , 

ipt + l'y^xxxx + bifixxx + ^J■Vxx + = 0, ip eC, {iy,b,fi) eTZ, ly ^0. (2) 

Sa réduction en onde propagative 

(p{x, t) = c + m(^), s, ~ X — et, (3) 

E{u, ^) EE vu'" + bu" + fiu' + y + A = 0, (j/, 6, /i) en, ly^O, (4) 

où A est une constante d'intégration, a également un comportement chaotique j21| . 
et elle dépend de deux paramètres sans dimension, /{pu) et vA/ . Nous noterons 
désormais x pour ^. Le problème à résoudre est le suivant. 

Problème. Trouver l'expression explicite de la "solution analytique générale" 
(définition supra) de l'EDO l@J. 

La motivation vient de ce que cette solution (pour CGL3 et donc pour KS) est 
"observée" sous forme de texture bien définie, tant dans des simulations numériques 
que dans de véritables expériences physiques [TT|. 

2 Suppression locale de la contribution chaotique 

Comptons d'abord le nombre (nécessairement inférieur à trois) de constantes arbi- 
traires dans la solution analytique générale. La recherche d'un comportement local 
algébrique au voisinage d'une singularité mobile xç, (mobile signifie : qui dépend des 
conditions initiales), 

u ~x^xo "oX^: wo 7^ 0, x = x-xo, (5) 
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conduit à la série de Laurent 



- I5bx-' + - ^) + 0(x°), (6) 



d'oià sont absentes deux des trois constantes arbitraires. Celles-ci apparaissent en 
perturbation 0], 



u = u 



(o)^e^(i)+£V2) + ..., (7) 




oii le petit paramètre e ne figure pas dans l'EDO Q. L'équation linéarisée en u^"^ 

est alors du type de Fuchs au voisinage de a; = xq, avec pour équation indicielle 
(q = —6 désigne le degré de singularité du premier membre _B de 

lim x^'^'^i^dl + uox'')x'^'' (9) 
= - 3)(i - A){j - 5) + 12{)v = v{j + - Vàj + 60) (10) 

-i^)(.-iî^)^o. (U) 

La représentation locale de la solution générale, 

u{xo,ec+,ec-) = 120i^x~^{Taylor(x) 

+£[c+x("+'^'/2Taylor(x) 

+ c_x(^^-*^)/2Taylor(x)]+0(e2)}, 

où "Taylor" désigne des séries convergentes de x, dépend bien de trois constantes 
arbitraires (xo,£c+,ec_) (l'indice de Fuchs —1 ne représente qu'une translation de 
Xo). Le branchement mobile dense provenant des deux indices irrationnels caractérise 
j31| le comportement chaotique, et l'unique moyen de l'éliminer est d'exiger ec+ = 
ec- = 0, c'est-à-dire e = 0, restreignant ainsi à un seul arbitraire la partie analytique 
de la solution. Le problème est donc de trouver une expression compacte pour la série 
de Laurent ©. 

3 Sur la suppression globale de la contribution chao- 
tique 

Par élimination de la constante mobile xq (donc de x— xq) entre la série de Laurent ((SJ 
et sa dérivée, la solution analytique générale satisfait à une EDO autonome d'ordre 
un, 

F(u',u) = Q, (12) 

que nous appellerons sous-équation de (0}, puisque cette dernière en est une consé- 
quence différentielle. 

Trouver la solution analytique générale équivaut alors à trouver une sous-équation 
autonome d'ordre un. En effet, si l'on admet que F est polynomiale et que sa solution 
générale est uniforme, alors 

1. (Hermite Vol. II, §139]) le genre de la courbe algébrique F = {) est zéro ou 
un, 
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2. (Briot et Bouquet pages 58-59]) sa forme nécessaire est 

m 2m — 2k 

F{u,u') = ^ ^ Œj^kU^u' =0, ao,m = l, (13) 

où m est un entier positif et les aj^k sont des constantes, 

3. (Briot et Bouquet) la solution est, pour le genre zéro, une fraction rationnelle 
de e"^ ou de x (a désignant une constante); pour le genre un, une fonction 
elliptique. 

Dans les deux cas, il existe un algorithme, implanté en Maple par Mark van Hoeij 
|12|. qui donne l'expression explicite de la solution. 

Malheureusement, à notre connaissance, on ne sait pas effectuer l'élimination de 
X — xo entre la série de Laurent et sa dérivée qui, si elle était possible, fournirait H12|l . 

4 Rappel sur la classe d'équations P{u^^\u) = 

P désignant un polynôme de deux variables à coefficients complexes, la classe d'équa- 
tions autonomes 

P(u("\u)=0, (14) 

indépendantes des dérivées d'ordre 2 à ti — 1 inclus, a été classiquement beaucoup 
étudiée. Soit W la classe de fonctions d'une variable complexe x 

W := {rationnera;), rationne^e"^), elliptique(a;)}, (15) 

011 a est une constante complexe. C'est en fait la classe des fonctions elliptiques et de 
leurs dégénérescences, car 

lim rationnel(e"^) = lim rationnel((e"^ — l)/a) — rationnera;). (16) 

a— *0 a— *0 

Les résultats connus sont les suivants. 

1. (Briot et Bouquet) Pour n — 1, si la solution est uniforme, alors elle est dans 
W. 

2. (Picard) Pour n — 2, toute solution méromorphe est dans W. 

3. (Eremenko [7]) Si n est impair, si le genre de la courbe algébrique P = est 
nul, et s'il existe une solution particulière méromorphe qui possède au moins un 
pôle, alors cette solution est dans W. 

Ce troisième cas (qui se démontre par la théorie de Nevanlinna) s'applique à 
pour b = ^ = 0, la solution est alors connue [H] 

6 = 0, ^ = : -60vp'{x - xo,0,g3), 52 = 0, 53 = , (17) 

lOoOi^^ 

p représentant la fonction elliptique de Weierstrass, et il en existe même une extrap- 
olation de codimension un ^H) 

b^^l6,.^0 : - = -60-P'-156p-g,.g. = ^,53 = ^^^. (18) 
Notre problème revient à extrapoler cette solution pour une valeur arbitraire de 
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5 Solutions de codimension non nulle 

Elles appartiennent toutes (du moins celles qui sont connues) à la classe W, Eq. H15|l . 
L'unique solution connue de codimension trois est rationnelle, 

5 = 0, /i = 0, A^Q l2Qv{x - ïq)"^- (19) 

On connaît six solutions de codimension deux toutes rationnelles en e^^ , 



2 32 X 19i^2 ■ 4 X 19j/ ' V 19 



— il — ôwk — tann — a: — xn) 

4 X 19i^/ 2 2^ ' 



f k k \^ f k k ^ ^ 

— 156 ( — tanh 77(2:^ — a;o) ) + 120;/ ( — tanh 77(2; — a^o) 

les valeurs permises figurant dans la Table Q 



(20) 



Table 1: Les six solutions connues de codimension deux de I^J. Elles ont toutes la 
forme (^0)1 . 



6V(a*^) 


vA/ 


vk'^ / ^ 





-4950/19^ 450/19-^ 


11/19, -1/19 


144/47 


-1800/47^ 


1/47 


256/73 


-4050/73^ 


1/73 


16 


-18, -8 


1, -1 



Enfin, l'unique solution connue de codimension un est elliptique, cf. (fH^ . et c'est 
une extrapolation de la dernière ligne de la Table 

Une belle propriété commune à toutes ces solutions est d'admettre la représenta- 
tion 



PLogcr, V^mv— + \^h— ■ ^ ^ 



dx ' 



(21) 



où a est une fonction entière dont l'EDO est facile à construire. 

Avant de rechercher plus avant la solution inconnue de codimension zéro, il con- 
vient de se poser la question si elle est ou non dans la classe VF, c.à.d. si elle est 
elliptique. 

On pourrait craindre en effet que cette classe W soit insuffisante, au vu d'une 
EDO qui ne diffère de Q que par le terme uu' P^ES] . 



E = -u" 9fcu" + (aiu - 26fc^)M' - 2aifcM^ ^ 2^k^u, 
et dont la solution sort de la classe VF, 
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U = — e p 
ai 



1 



Xq, 92,93 , (^0,52,53) arbitraires. 



(22) 



(23) 



Une solution dans cette nouvelle classe est exclue pour Ql, car l'EDO d'ordre un 
résultant de l'élimination de entre (|23|l et sa dérivée, 

(u' - 2kuf - 4^3 + e'^^'-'92U + e^'^'^gg = 0, (24) 

n'est autonome que pour la dégénérescence rationnelle 92 = 93 = 0. 



6 



6 Vers une approche numérique de la solution el- 
liptique 

À partir de l'unique élément d'information que constitue la série de Laurent ((HJ, nous 
allons définir plusieurs approches numériques pour tester la validité de l'hypothèse 
elliptique (appartenance à la classe W). L'outil naturel est celui des approximants de 
Padé (AP) El- 

Etant donné les N premiers termes d'une série de Taylor en a; = 0, 

N 

5w = ^c,a;^ (25) 

j=o 

l'approximant de Padé [L, M] de la série est défini comme l'unique fraction rationnelle 

[i>^^] = 75K^>&o = l, (26) 

satisfaisant à la condition 

Sn - [L, M] = 0{x^+^), L + M = N. (27) 

L'extension aux séries de Laurent ne présente pas de difficulté. En particulier, les AP 
sont exacts sur les fractions rationnelles dès que L et M sont suffisamment grands. 

Le calcul préalable d'un grand nombre (environ 200) de termes de la série de Lau- 
rent se fait par un algorithme linéaire, la seule difficulté, aisée à surmonter, étant de ne 
jamais engendrer de termes qui s'avéreront inutiles. Il importe également de calculer 
les Cj, puis les (a/, 6m), sur Q, jamais sur TZ, à cause d'importantes compensations. 
C'est pour cela que nous avons dû écarter le programme pade de Numerical recipes 
P7| . qui engendre des doublets de Froissart [21 El indésirables (couples d'un zéro et 
d'un pôle très proches, ces couples se répartissant sur un cercle sans définir de coupure 
comme ce serait le cas s'il ne s'agissait pas d'erreurs numériques). La fonction pade 
du langage formel Maple est en revanche tout à fait adaptée. 

Remarque. Il est possible d'éviter l'étape de calcul des AP et d'obtenir directement 
et rapidement PO] les valeurs numériques des zéros et des pôles du Padé [L, M]. 

6.1 Premier test d'ellipticité 

La structure des pôles et des zéros des AP de © caractérise la classe de fonctions à 
laquelle appartient sa somme inconnue. 

Une somme elliptique est caractérisée par une double infinité de pôles et de zéros 
situés aux noeuds d'un réseau doublement périodique, une somme rationnelle en e"^, 
avec a complexe, par une simple infinité de pôles et de zéros, une somme rationnelle 
en X par un nombre fini de pôles et de zéros, etc. C'est pour tester les structures 
doublement périodiques que l'on a besoin d'un grand nombre de termes dans la série 
de Laurent. 

Pour tous les cas de solution connue f (|20|l . I|18|l ). en choisissant pour les constantes 
fixes (z^, 6, /i,74) des valeurs rationnelles, le résultat est conforme aux espérances, voir 
Fig- n Pour plusieurs cas génériques de valeurs de b^/ipiy) (par exemple b^/{^iv) = 
4, 8, Fig.lSJ, nous trouvons une solide indication numérique que les textures observées 
sont toutes doublement périodiques, ce qui est la signature des fonctions elliptiques. 
Dans tous les cas, connus et inconnus, il est aisé de compter le nombre (commun) de 
zéros et de pôles par période, et les calculs indiquent que ce nombre serait égal à trois. 

Remarque. La connaissance du nombre et des positions des zéros et des pôles 
d'une période permet la détermination numérique des invariants g2 et g^. 
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Padé (Laurent e Q) [30,30] 



. (a) 



Exact (Half Periods = ( os,, os; ) ) 



. (b) 



half period 



.?<" O .X-' O 



b= 16 
11=15 
A = -1408 

S, = 48 



b= 16 
H= 1& 
A = -]40K 



) p' - 240 jJ - 64, = 64/3. = 



Errors: Padé (Laurent e fî) [100,100] 



Padé(Laurent e ()) [30,30] 



_(d) 



^l- 16 

A = -1408 
g,. 64/3 
g, = 48 

A<0 



b= 16 
11= 16 
A = -52168 



Figure 1: Plan complexe x du seul cas elliptique connu (fe^ — A = arb). 

La structure des zéros (symbole 0) et des pôles (symbole X) doit être doublement 
périodique, avec trois zéros simples et un pôle triple par parallélogramme période. 
Figure (d) : {v,b,n,A) = (1,16,16,-52168), calcul de l'AP [30,30] sur Q. Les 
trois autres figures sont tracées pour {i/,b,fi,A) = (1,16,16,-1408), donc (52,^3) = 
(64/3,48). Figure (b) : emplacements exacts des zéros et des pôles, calculés à partir 
de la solution exacte lfTH|) . Figure (c) : zéros et pôles de l'AP [100,100] calculé sur 
TZ, montrant l'accumulation des doublets de Froissart sur le cercle unité (échelles 
différentes sur les deux axes). Figure (a) : zéros et pôles de l'AP [30,30] calculé sur 
Q. Dans tous les cas, on discerne nettement les trois zéros simples et le pôle triple de 
chaque période. 



6.2 Deuxième test d'ellipticité 

Si la solution u est elliptique, alors, d'après un résultat classique, il existe deux frac- 
tions rationnelles Ri,R2 telles que u = Ri{p) + ^2(p) p' , plus précisément 

^Polyg(p)^ Poly^(p) 
partie paire 

oit Poly^) désigne un polynôme de degré D. 

En vue d'une recherche directe de la solution sous la forme (|28|l . il est en principe 
possible de calculer numériquement les degrés Ni,N2,D des fractions rationnelles de 
p, en remplaçant le plan complexe de x par celui de p{x), de la manière suivante. 

Tout d'abord, la nécessité d'une seule variable complexe dans les AP oblige à 
éliminer p', donc à introduire le branchement algébrique indésirable Y^4p'^ — 52 p — ^3 ■ 
Un moyen simple de le supprimer est alors, puisque p est une fonction paire, de se 
restreindre à la partie paire de la série de Laurent u{x = x — xq), 

^ ^ 19 X 321.2 192 X 64Z.4 ^ + ^' ^^^> 

d'inverser vol. II, p. 527] la série de Laurent p de x, 

pix, 92,93) = + ^x' + + 0(x«) (30) 
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(a) 



Padé(Laurent ^ Q) [20,20] 

"1 ' \ ' T ' \ ' \ ' 

\ O 

O Q o é "9 

M o ^ o ^ o 
,6 d o .o \. 



o 














= 1 


b 


= 1 




= 1/8 


A 


= 1000 



(b) 



Padé (Laurent e Q) [20,20] 



- o o 



(c) 



Padé (Laurent e g) [20,20] 



: o é d O 
^ ■■■îj' o x.' o .-(x 

o 'o .o "P, .é 



■■■""o"- . 



o 



2 4 



A = 1000 



Padé(Laurent e g) [30,30] 



(d) 



^ ' I o ^ 

>o o o 6 o ® ,, 
>o ° ,6 P .o ® 

* o Si; o 



-2 



11.1/4 
A = 1000 



Figure 2: Plan complexe a; pour les valeurs h"^ /{jiv) = 4 et 8 (cas génériques inconnus) 
et desAP calculés sur Q. Fig. (a), {u,b,^i) = (1, 1, 1/8), v4 = 1000. Fig. (b), {v,h,^i) = 
(1, 1,1/8), A = -1000. Figs. (c) et (d), {v,h,ii,A) = (1,1,1/4,1000), pour [L,M] = 
[20,20] (Fig. (c)) et [30,30] (Fig. (d)). La structure est, avec une bonne précision, 
doublement périodique à trois zéros simples et un pôle triple par période. 



en la série de Laurent de p. 



.92 2 

+ 20^ 

( 9i 



55i 



75i ,-4 



P 



293233 



+ 



V 12480 1274 

4 



P 



6693253 



1200 " 3080 
; 167g|g3 , 
73920 ^ 

r« + 0(p-9) 



^509184 340340 
puis de composer les deux séries (|29|l et (|31|l pour obtenir 

6(56/11/ - 135^) / 36.92 b {iQiiv - 3P 



(31) 



paire 



-156p 



19 X 32j/2 



192 X 64;/4 



P 



(32) 

La question est alors de tester si cette série H32|l . dont les coefficients dépendent 
de {v, b, /z. A, (72, 93), a bien pour somme une fraction rationnelle 



paire 



(a;, 52, 53) = 



Polywi (p) 
PolyD(p) ■ 



(33) 



La difficulté provient de la nécessité de donner, dans la série des valeurs 

numériques, de préférence rationnelles, à {92,93) avant de pouvoir calculer les zéros 
et les pôles des AP. Puisque le parallélogramme période du plan x est déjà déterminé, 
il permet de calculer des valeurs approchées (complexes) de (32,53). H importe donc 
d'étudier la sensibilité du résultat (valeurs de Ni,D) aux valeurs de (52,53), et cette 
sensibilité ne peut être testée que sur l'unique solution elliptique connue H18(l . La 
Figure 13 montre la structure des AP de dans le plan complexe p pour 6^ — lôfii' 
et pour un choix de (52, 53) voisin de celui de ce cas. 
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Padé (Laurent e Q) [L,M] 



Padé (Laurent e Q) [10,10] 



(a) 



(b) 



b = 16 

|l= 16 
A = - 140S 
g, = 64/î 
^48 



v = l 

p. 16 
' A = -14(18 

E,-45 



u - - 15 b - / (4v) 



Figure 3: Plan complexe p de la partie paire de la solution u pour b'^ /(fiv) — 16, A = 
— 1408, ce qui définit par l(TÏÏ)l les valeurs 32 = 64/3,^3 = 48. Dans la figure (a), 
{iy,b,fj.,A,g2,g3) = (1,16,16,-1408,64/3,48), tout AP [L,M] est exact car calculé 
sur Q, la figure contient donc l'unique zéro p = —fi/ifiOu) et elle n'a aucun pôle. 
Dans la figure (b), où les valeurs de (52,53) sont choisies rationnelles et légèrement 
différentes des valeurs exactes, (1^,6,^,^,52,53) — (1,16,16,-1408,22,49), l'AP 
[10, 10], calculé sur Q, fait constater 10 quasi-compensations d'un pôle par un zéro. 



6.3 Détermination des degrés de la sous-équation 

Suivant une suggestion de Belabas faite lors de la conférence, les degrés de la courbe 
algébrique 1)13(1 pourraient se calculer par une technique de Padé, comme suit. A partir 
de la série de Laurent on construit d'abord une autre série de Laurent à pôle sim- 
ple, par exemple u' /u — Laurent (a; — a;o), on l'inverse enx — xo — Laurent (m' /u), et on 
la substitue dans la série u — Laurent(a; — a;o). La série résultante u = Laurent(it'/u), 
qui réalise sous forme locale l'élimination de x — xq entre u et u' , est alors analysée 
par Padé. Du branchement algébrique ainsi trouvé, il devrait être possible de déduire 
le degré en U de la courbe algébrique G{U, V) — 0, avec U — u,V — u'/u, donc les 
degrés de lfTÏÏ)l . 



7 Approches analytiques 

L'étude numérique précédente semble conclure à la nature elliptique de la solution 
cherchée. Nous présentons maintenant brièvement trois méthodes analytiques sus- 
ceptibles de fournir les coefficients de la sous-équation H12|) . Si de plus celle-ci est 
supposée algébrique et à solution générale uniforme, alors elle a la forme nécessaire 

m- 

7.1 Équation de Monge- Ampère équivalente 

La fonction inconnue de deux variables F(u, u') qui définit la sous-équation d'ordre 
un ((12|l satisfait à une équation aux dérivées partielles très simple, obtenue par 
élimination de u" et u'" entre et les deux équations 

^F{u{x),u'ix))=0, (34) 

j2 

-^F{u{x),u'{x))^0. (35) 



Dans la notation classique 



dF dF 

Q-ip d'^F Q-2p 

^ = â^' ^ = dxdY' ^ = âï^' ^^^^ 
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c'est une équation de Monge- Ampère qui s'écrit 

Y^i-Q^R + 2PQS- P-T) + YP^Q _ ,yP ^ ^ ^ ^ ^ (3,) 

Notre problème équivaut alors à en trouver une solution particulière à zéro paramètre 
qui soit polynomiale. Ce nouveau problème pourrait n'être pas plus simple que 
l'original, car la méthode des caractéristiques est inapplicable puisque l'équation (|37|l 
hérite de la nature chaotique de KS. En particulier nous n'avons pas réussi à obtenir 
les degrés de cette solution polynomiale, si elle existe. 



7.2 Construction de la sous-équation par l'algorithme LLL 

Il existe une méthode [201 héritée de la cryptographie, qui fournit en nombres 
entiers les coefficients de pour m donné (et ao,m non normalisé à un). Connue 
sous le nom algorithme de réduction du réseau (lattice réduction algorithm) ou 
algorithme LLL, elle reçoit en entrée l'entier m, la suite des inconnues aj^k de H13|l . 
et L valeurs numériques {ui,u'j),i = 1, • • • ,L (donc dans TZ et non pas dans Q) de 
la somme des séries de Laurent u et u' aux points Xi — xo,i = 1, ■ • • , L. En sortie, 
si L est suffisamment grand devant le nombre (m + 1)^ d'inconnues aj,k, elle exhibe 
en un temps polynomial une suite d'entiers aj,fe, par résolution du système linéaire 
homogène surdéterminé 

rn 2m — 2k 

V/ = 1,...,L : J2 E aj,kuiu['' ^0. (38) 

k=0 j=0 

Si les entiers ainsi trouvés sont tous nuls, il faut augmenter m et recommencer. Nous 
n'avons pas encore mis cette méthode en œuvre car celle de la section suivante semble 
beaucoup plus prometteuse. 



7.3 Construction de la sous-équation par la série de Laurent 

Deux d'entre nous ont récemment proposé j23| une méthode analytique permettant 
de construire la sous-équation H13(l à partir de la série de Laurent L'algorithme 
est le suivant. 

1. Choisir un entier positif m et définir l'EDO de Briot and Bouquet 1)13(1 . Elle 
contient (m + 1)^ — 1 constantes inconnues Oj^k- 

2. Calculer J termes de la série de Laurent, où J est légèrement supérieur à (m + 
1)^-1, 



avec dans le cas de KS p = — 3. 

3. Exiger que la série de Laurent satisfasse à l'EDO de Briot and Bouquet, c.à.d. de- 
mander l'annulation identique de la série de Laurent de F{u, u') à l'ordre J 



F = ^J^' + j , D^mip- 1), Vj : F, = 0. (40) 

S'il n'existe pas de solution pour aj.fe, augmenter m et revenir à l'étape 1. 
4. Pour toute solution, intégrer l'EDO autonome d'ordre un (fT^ . 
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Le cœur de la méthode est la troisième étape, où le système de J + 1 équations 
Fj — aux (m + 1)^ — 1 inconnues Oj.fc est linéaire and surdéterminé, donc très facile 
à résoudre. 

Quant à la quatrième étape, elle est résolue par le progiciel algcurves 12 , qui met 
en oeuvre des algorithmes de Poincaré. 

Dans le cas de l'équation Q étudiée ici, du fait de la valeur p ~ —3, la valeur 
minimale de m est m = 3, ce qui correspond à 15 inconnues Uj^k et même à 10 
inconnues seulement, à cause de la règle de sélection supplémentaire m{p — 1) < 
jp + k(p — 1), qui provient du degré de singularité p = —3 de u. Ces dix coefhcients 
sont déterminés par le système de Cramer défini par les dix équations Fj = 0, j = 
0:6,8, 9, 12. Le système résiduel, qui est non-linéaire et surdéterminé aux inconnues 
fixes (i/, 6, /i, A), a pour PGCD — l&jiv (cf. lfTÏÏ)l l ce qui définit la solution de 
codimension un, 



— = 16, U' + — u J lu'-—uA+—[ul + ^ + — ] = 0, 
fil' \ 2v } \ 4:1/ J 40i/ \ 1/ S J 

36/i 

Us ^ U + — 

2v 



(41) 



Après division par ce PGCD, le système restant aux inconnues {i>,b,fj,,A) admet 
exactement quatre solutions (il est suffisant d'arrêter la série h j — 16 pour obtenir 
ce résultat), 

6 = 0, 

180 Y / , 360Ai2\ 9 / 2 30/^ , 30V^\^ „ 

,3 9 / 2 30Ai , 30VA^ 
b = ^,u'' + —U + -^u' + -£-] =0, (43) 



AOv \ 19 193i^ , 

144 56^ / , 6 \^ 9 4 
= u :— , U + — u., + — — Uc 



fiv 47 47v \ Av J mv 

6^ _ 256 _ 456Ai 
^ ~ 73"' 584^' 



= 0, (44) 



2 



b \ ( , b \ 9/2 56^ 56^ , , 

H Us u H Us H ui H -Us H u' = 0. (45) 

8v ) V 2zy 7 40i/ V 10241/2 128i/ / ^ ' 

A l'étape 4, on trouve alors la valeur un pour le genre de (|1T)| . et la valeur zéro 
pour celui de puis les expressions fTÏÏ|l et (^0)1 pour les solutions. 

Avec la valeur minimale m = 3, on retrouve donc tous les résultats connus. Le 
calcul pour m > 4 est actuellement en cours. 



8 Conclusion 

L'étude par approximants de Padé semble conforter la nature probablement ellip- 
tique de l'unique solution analytique de codimension nulle de l'EDO chaotique 
Parmi les autres pistes, mi-analytiques, mi-numériques, qu'il serait possible de met- 
tre en oeuvre, il en existe une particulièrement intéressante car constructive, c'est 
celle d'obtention de la sous-équation par la série de Laurent. Son seul inconvénient 
est l'absence de borne supérieure pour le degré m, ce qui contraint à de volumineux 
calculs qui sont encore en cours. 
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